clear all;clc

%% 简单线性回归
x = 1:100; 
y = -0.3*x + 2*randn(1,100); 
[p,S] = polyfit(x,y,1); 
f = polyval(p,x); 
%% 误差估计，其实在预测上并没有任何区别！！
[y_fit,delta] = polyval(p,x,S);
plot(x,y,'o',x,f,'+',x,y_fit,'*') 